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Abstract. We approximate the spectral data (eigenvalues and eigenfunctions) of compact 
Riemannian manifold by the spectral data of a sequence of (computable) discrete Laplace 
operators associated to some graphs immersed in the manifold. We give an upper bound on 
the error that depends on upper bounds on the diameter and the sectional curvature and on 
a lower bound on the injectivity radius. 



1. Introduction 

We prove that the spectral data (eigenvalues and eigenfunctions) of any closed Riemannian 
manifold can be approximated by the corresponding spectral data of the Laplace operator of 
some graphs geodesically immersed in the manifold. It is an extension of the finite elements 
method to the Riemannian setting. The two main points of our method are the following. 

(1) The error made on the spectral data are bounded above by universal functions of some 
bounds on the geometry of the manifold (i.e. bounds on the injectivity radius, the 
sectional curvature and the diameter) and of the graph (i.e. bounds on the thinness 
and mesh of the graph). This errors tend to as the mesh of the graph tends to 0. 

(2) The Laplace operator of a metric graph is a universal and explicitly computable function 
of its edge- lengths. 

Before stating our main results, we need a few definitions and notations. 

1.1. Definitions and notations. We will work with a special kind of immersed graphs, that 
we will call geodesic triangulations (see the definition in Section They are not necessarily 
actual triangulations of M (for instance the simplices of dimension greater than 1 are not 
necessarily realized as subset of M) but are more easier to construct. 

A geodesic triangulation T of an n-dimensional Riemannian manifold [M'^^g] is a set of 
points {xi)i^i^iq of M endowed with a structure of abstract simplicial complex K which satisfies 
the two properties 12.11 and 12 . 21 of section [21 We denote by Sp the set of closed p-simplices of K. 
We identify the edges of K with some minimizing, geodesic segment between their vertexes. 
For any closed simplex a of K, we set St{a) (resp. Stp{a)) the set of the closed simphces (resp. 
of dimension p) of K that contains a. The vertices of any a € Sp are naturally ordered by 
their indices {a — {0:^^(0), . . . , Xj^(p)} with v(0) < . . . < ia-ip))- We set X^- ~ •'2^j„(o) and for any 
distinct < fc < n, we set a vector of Tx„M such that Xi^{k) — exp;,^-^ (w^). We set also A'^ 

the associated Gramm matrix ( g{v1, Vi) ) 1 < ; < „ • Given a geodesic triangulation of (Af", g), 

1 < fc < n 

we note rriT its mesh (the maximal length of its edges) and Ot its thinness, i.e the quantity 

^ / , , length(ei) 

Bt = max max mTldetA ) , max — — - 

\ a e Sn {ei.e2)eSi length(e2j 
< fc < n 



Key words and phrases. Spectral theory, graphs, finite elements, Riemannian geometry, finite elements, 
discretization. 
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Eventually, on the set of functions j/ : T — > R, (where N is the number of vertices of T and 
we identify T with Sq), we define two quadratic forms by the formulae 

(1-1) 1^'^ = ?;^ ^ y-y^ ^ 

1 " 

(1.2) qriy) ^^ Y. ^cletA- ^ (A")'''(y*.(fc)-2/«„(o))(y»„(0-y..(o))- 

Note that if is a sub-complex of M" then \y\T and (7t(2/) give respectively the L^-norm and 
Dirichlet energy of the affine-by-parts expansion of y. 

1.2. Main results. For any closed, Riemannian n-manifold {M,g), we denote by Sm its di- 
ameter, by R an upper bound of all its sectional curvatures and by iM its injectivity radius. 
We denote also by = Ao(T) < • • • < AAr_i(r) the eigenvalues of qr with respect to | • |y and 
= Xo{M) < \i{M) < ■■ ■ < XiiAI) < ■■■ the eigenvalues of {M"',g). 

Theorem 1.1. Let n > 2 be an integer, and e g]0, l[ be a real number. There exists a constant 
C{n) such that if 

i) {M",g) is a closed, Riemannian n-manifold which satisfies \R\ < A^, 

ii) T is a geodesic triangulation of AI which satisfies < C{n){ ^ o^'ec^p ) ^' 
then we have 

(1 - e)Ap(r) < Ap(M) < (1 + e)Ap(T). 

Remark 1.2. The constant C{n) is computable. The second condition says that any finite 
number of eigenvalues can be approximated provided the mesh of the graph is .small enough and 
the thinness is controlled. 

Remark 1.3. The matrices A"' depend on the angle between some edges of T issued from a 
same vertex, but the same result is valid if we replace the coefficient g{v^,vf) by 

i d'^{X„,x,^(^k)) + d^{X„,x,^^i)) - d^{x,^(^k),x,^(,^) 

in the definition of the matrix A'^ . This gives approximation of the eigenvalues of M by the 
eigenvalues of a discrete Laplace operator whose coefficients are universal functions of the 
lengths of a geodesically immersed graph of M . 

Note that in [4 , the authors get the same result for another geometric quadratic form qx, 
whose coefficients depend on the volume of the Voronoi' cells associated to a lattice {xi)i^i which 
need not to be the vertices of a geodesic triangulation. 

We denote by {ff) the eigenvectors of qr with respect to | • and let {fi)i£t>i be a 
orthonormal family of eigenfunctions of {M'^,g) such that A/i = A^/i for all z G N. For 
some integers p < q, we set Ep^g (resp. Fp,q) the sum of the eigenspaces of A(M) (resp. qr) 
associated to the eigenvalues (Ai(M))p+i^i^q (resp. (Ai(T))p+i<ji^^ and Pp^q (resp. Qp^q) the 
normal projection on Ep^q (resp. Fp^q). 



Theorem 1.4. Under the assumptions of Theorem if there exist some integers p < q 
and rj > such that Xp + rj ^ Ap+i and Xq + ri ^ Aq+i, then for any f £ Ep, we have 

C(o n A ) 1 

||i?(/)-Pp,,oi?(/)||2, sc: (f^)^||fl(/)||j-, and for any (y,) G Fp, we have \\W{y,)~ 

' rriT 



]p,, o iy(y,)||| C{q,n,A, ,^){f^)^ \\W{y,m 
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To get an approximation of the spectral data of (M", g) by those of qr in Theorems 11.11 and 
11.41 we need some geodesic triangulations with arbitrary small mesh but bounded thinness. 
The existence of such fat triangulation is often admitted or conjectured but we do not know 
complete published proof of this fact. For sake of completeness, we give a constructuve proof 
of the following result (based on the previous work of J.Cheeger, S.Miiller and R.Schrader). 

Theorem 1.5. Let n > 2 be an integer and D^iq and A be some positive real numbers. There 
exist some constants l3{n) and C{io/D,A,n) such that for any Riemannian manifold {M'^,g) 
with diameter Sm ^ D, sectional curvature Sm \R\ ^ and injectivity radius iM ^ io and for 
any e s]0, C(io/-D, A, n)[, there exists a geodesic triangulation T of M with mesh mr < e and 
thinness Qt < l/C(io/I?, A, n). 

Remark 1.6. Once again the constants of Theorem ] 1.5\ are explicitly computable. Combining 
Theorems ] 1 . 5\ and \l.l[ for any compact manifold {M"',g), any iV g N and any e > we get a 
method to construct a geodesic triangulation T of M such that we have (1 — e)Ap(T) < Ap(M) < 
{l + e)Xp{T) for any p ^ N . 

1.3. Main steps of the proof. Let (i?, (• , •)) be a Euclidean space endowed with a bilinear 
symmetric form q, and Ao < • • • < XdimE-i be the eigenvalues of q with respect to (•,•). Using 
the min-max principle we readily infer the following spectral comparison principle. 

Proposition 1.7 (small eigenvalue principle). Let (•, and (i?2; (•, •)2) be two Euclidean 
spaces endowed respectively with quadratic forms qi, 52 • If there exists a linear map ^ : Ei ^ E2 
and two positive real numbers a, /3 such that 

{(^{x),(^{x))2 > a{x,x)i and q2{^{x),^{x)) < P qi{x,x), 

then we have Xk{q2) < a^kiqi) for any k. 

Proposition II . 71 is usually used to compare spectra under small perturbations on the metric 
or on the manifold. It is the key tool of our eigenvalues approximation method. 

Given a manifold M and a geodesic triangulation T of M, we denote by {xi)i<i<N the 
vertices of a T, by Ep the subspace of H^'^{M) spanned by the p + 1 first eigenfunctions 
{fi)o<i<'p of M , hy R : Ep ^ the natural restriction map R{f) = {f{xi))i<i<N , by (•, •) the 
scalar product on Ep induced by the L^-norm on M, and we set q{f) = Jj^j \df\'^. The spectrum 
of q with respect to (•, •) is given by the p + 1 first eigenvalues (counted with multiplicities) of 
M. We then proceed in two steps. 

(1) A Moser's iteration scheme gives bounds of the quotients on Ep \ {0} by a 

universal function of Ap, S and A (see proposition 13. ip . This Hessian bounds imply the 
following estimates (see Propositions 14.11 and 14. 2p 



(1-3) {R{f),R{f))T- / / 



M 



< Cmr / 

M 



(1.4) 9T(i?(/)) < {l + CmT} / \df\\ 

for any f G Ep, where | • |t and qx are the discrete quadratic forms on given by the 
formulae (jl.ip and (jl.2p . and where C is a constant which depends on bounds on Ap 
and on the geometries of (M", g) and T. Proposition 1 1 . 71 gives then some lower bounds 
on the spectrum of (M",^) of the form (see Theoremg^J) Afe(M) > (l - CmT)\k{T) 
for any k < p. 
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(2) In Section 231 we construct an expansion (Withney) map W : — > C°°(Af) such 
that RoW = Wrn and which satifies the following estimates 

(1-5) UvMt- f W{y,f <CmT{\\{y^)fT + qT{y^)), 

(1.6) / \dWiy,)\^ <{l + CmT)qM, 

for any (j/i) e R^. From Proposition 11.71 again we infer that for any k < p we have 
(l — CniT) Xk{M) < Afc(T) (see Theorem 14. 8p . The construction of the withney map is 
the main technical difficulty of the proof. It is done by local mean of the affine expan- 
sions obtained by identifying the simplex of the geodesic triangulation with Euclidean 
simplicex through the Riemannian exponential maps at the vertices of the simplex. 

Note that J. Dodziuk ^ developed another generalization of the finite element method to 
compact Riemannian manifolds in which, to any smooth triangulation of (A/, g) is associated 
the subspace of i7^'^(M", g) of the continuous functions on M which are affine on each simplex 
(this subspace has finite dimension). The authors consider on it the quadratic forms induced 
by II • II 2 and the ambient Dirichlet form q{f) = Jj^^ Id/p . They prove that the spectrum of the 
discrete Dirichlet form with respect to the the discrete norm converges to the spectrum of 
(M", g) when the mesh of the triangulations tend to with controlled thinness. However, they 
do not prove that the error is bounded by geometrical bounds on the manifolds, and moreover, 
the discrete quadratic forms cannot be explicitly computed as function of the geometric data 
(edge's lengths, edge's angle) of the triangulations. 

The proof of Theorem 11.41 is done in section [5] using the above estimates and the technique 
developed by Y. Colin de Verdiere in [5] . 

Aknowledgement We thank S.Gallot for fruitful discussions and C.Vernicos for bringing our 
attention to the paper [1]. 



2. Geodesic triangulations 

2.1. Definition. A geodesic triangulation T of an n-dimensional Riemannian manifold {M"-,g) 
is a set of points (xi)i^i^Ar of M endowed with a structure of abstract n-dimensional simplicial 
complex K whose simplices are all contained in a n-dimensional simplex of K and which satisfies 
two more properties for which we need to complete the notations of the introduction. 

For a G Sp and any distinct < k,l < p, v^i is a vector of Tx-^^f.^M such that a;i^(;) = 

exP:r.„(.)K/) and \v'^i\ = = d{x,^(^k), x^^(i)). We set the matrix (g(w^,,w^„)) ^ ^ ■ 

m ^ k 

For any vertex a; of T and any a £ St{x), Ca is the cone of TkM spanned by the vectors 
i'"lf)-^(x)i)o<i<P- l^ote that for any a € Sn, we have = a::j^(o), = L'^^. and A" = A^ 
(according to the definitions given in the introduction section). For any simplex a of if, we 
denote by Np{a) the set of all the simplices of dimension p that intersect a. 

We set F = {{9,) e R"+7 Y^^0i = 1}. For any A G M, is the closed n-simplex F f] 
[1 — A, -|-cx)["+^ (we will denote A" = A"). Given a G Sn and < k < n, we get some local 
barycentric coordinates on M by the formula 

B^,:{9i)eF>^exp,^^^^^(j2Siv^ki)- 

l^k 

We set A • To- = -Bq (A"). Eventually, a geodesic triangulation T has to satisfy the following 
two properties. 



(2.1) For any vertex x of T, {Ca)a-est(x) induces a triangulation of the unit sphere of T^M. 
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(2.2) For any disjoint ct, a' e 5„ and any < A:, fc' < n, we have Bl (A") n B'(, (A") = 0. 

Eventually, a geodesic triangulation with boundary T of an n-dimensional Riemannian 
manifold (M^^g) is a set of points (xi)i^i^Ar of M endowed with a structure of abstract n- 
dimcnsional simplicial complex K whose simplices are all contained in a rt-dimensional simplex 
of K and which satisfies condition (|2.2p but condition (j2.ip only for vertexes not on the bound- 
ary of where we call boundary of K the complex of the simplices of K that are contained 
in a n — 1 simplex of K itself contained in only one n-dimensional simplex of K. 



2.2. Metric estimates. We now study some metric properties of the geodesic triangulations 
in bounded geometry. We first recall some estimates on the Riemannian exponential map whose 
proofs can be found in [|^5,. 

Theorem 2.1. Let {M",g) be a compact, Riemannian manifold with (5|,/|cr| < A^. 

Let V e T^M be fixed, y = exp^(u) and for any w € T^M , let w(t) be the parallel translation 
of w along 1 1— >■ exp^(to). // we define two maps from T^M to M by F{w) — exp^(-y + w) and 
G{w) — expj^ , then they satisfy the following estimates 

d{F{w),G{w)) < ldix,y)\w\g^^smh{^idix,y) + \w\gj), 
6 Om Om 

sinh(AM) 

K exp.H - w{l)\ < I^K— - 1). 

Sm 

Theorem 2.2. Let (A/", 5) be a manifold with S1j\(j\ < A^ and e < mi(iM, ^) be a positive 
real. Then for any x G M , the map exp^, is a diffeomorphism from B{Ox, e) C T^AI to B{x, e) 
and for any u,v Cz B{Ox, e) we have that 

{l-A\^f)\u-v\,^<d{c^p^{u),exp,iv)) < {l + A^{^f)\u-v\,^. 
OM Om 



If r is a geodesic triangulation of (Af",f/) with mesh smaller than ii\i /lO, then for any 
a G Sn and any < k < n, the map gives some coordinates on a neighbourhood of A"q. 
We can compare these coordinates for fixed cr but different values of k. 

Lemma 2.3. Let {M",g) be a manifold with Sf,j\a\ < A^ and T be a geodesic triangulation 
with mesh lOmx < inf (im, ■ '^^V ^ '^"■2^ < fci, ^2 < n and any {9i) g A"o, we 
have 

d{Bl{e,),Bli9,)) < lOA^'^fmrY^m. 

Proof Let {9,) e A^q. We set v = vl^^^ and we = E; Oi{vl^i - Then we have = 

expj.^ ~^ "^eit) is the parallel transport of w along s 1— >■ exp^,. ^t,^^{sv%_^^^) then 

Theorem 12.11 implies that 

d(B^^(0O,exp,^^^^^^(u;e(l))) -d(exp,^^^^^^(i; + w;0),exp,^^^^^^(u;e(l))) < {-^frn%Y,\e,\. 

i 

For = 5ii we get 

^K.('5.0,exp,^^,^^,(w;5jl))) -d(exp,^^^^^,(z;^^;),exp,^^,^^,KJl))) < 
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SO, by Theorem l2.21 we have \u'^^i — v'^^i ^ < 2{j^)'^m^, where u^^j is the parallel transport 
from Xi^^ki) to x,^(^k2) of the vector w^^, - w^^j.^. Hence we get 

d{B^,^ie,),Bli9,)) < d(S,-^(^?,),exp,^^,^^^(we(l))) +d(exp,^^^^^^(we(l)), 5^,(0.)) 



OM Om , 



□ 



For any t Sp, with l<p<n— l,we set 

n = U exp..^,,, ({E > 0, E ^ 1} 



creSt„(T) j=o 



T.= U ULoexp.^^,,,({E^.<.-(..^,,)/^.>0,E^.^l})- 

creSt„(T) J=0 

For any subset A C M, we set B{A, r) the tubular neighbourhood of A and radius r. If ^ is 
empty we set B{A, r) = 0. 

Lemma 12.31 implies the following result. 

Corollary 2.4. There exists a constant C{n) such that if [M" , g) and T satisfy 5\i\a\ < 
and lOmx < inf (zm, c(n')e^"A ) ^^s*^ have the following properties 

(1) St^ < mT(l + 10(3^)2m|) /or any a G 

(2) i/ie ((1 — jy) • 117)^^5 dfe disjoints and the ((1 + 77) • T^<j) cover M , 

(3) /or 077.?/ a,T € K such that a H t = $, we have d(T^, T^) > ^("^93^; 

(4) for any a € Sp and t G K, the tubular neighbourhoods _B^Tr, o^^^-^"^T^ and B[Ta\ 
B{T^nT,aPf^mT),aP+^f^mT) are disjoint, 



where 7] = C{n) B|?^(^)^to|. aTid a — (ji^^^q-i^ ■ 

Proof. We set 7Ti = 7TiT and O = Qt- Since the points Xi^(^k) are in B{X^,m) , the n-simplex 
Aa = (Ox„,woi' ■ • ■ '«OrJ of Tx^M has a diameter less than ^_fJ^-.2^2 < w(l + 2{-^fm'^). 
Moreover, if a' is a (ri — l)-face and H is the iso-barycentre of Ao-, then the distance from H 
to cr' is equal to („"^^°\^fq.' — c(J)e^" ■ Hence, Theorems 12.31 and 12.21 imply that wc can choose 
C{n) large enough so that T^r have diameter less than 777(1 + 4(^)^777^), the (1 + 77)Tct contain 
all the B^(A") (0 < fc < 71) and (l-77)T^nE^'(A") = for ah a' e 5„\{ct} and ah < fc < 77. 
In Particular the ((1 — v)Tcr)^^g are disjoint. 

Let T e S'p be a simplex of T. We now show by recurrence on p that 

Nr^ 11 B''^^, ,(A") 

a e N„{t) 
<k <p 

is a neighbourhood of Tt- in M and that d{Tr,dNr) > (g^„™Q^„ . 

Note that for any a £ Sn, A^- has heights greater than > ^t^;^^^^, where a' is the 

face of a with dimension n — 1 and smallest volume. So the case p = derives from the first 
axiom of geodesic triangulations and from Lemma 12.31 
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If T is of dimension p> 1 then exp^,. {{^2^=0 )/ > 0, ^ 6j < 1}) is interior 

to (J ^^-^{i (J,-,-) (A") for any < fc < p, by the first axiom of the geodesic triangulation 

o-e5t„(3;i^(fc)) 

and the smallness of m. Its boundary is also included in by the recurrence assumption and 
Lemma [2.31 So is a subset of Nr- It remains to show that d{T^,dNT) > p^„™g)2n , which 
combined with Lemma [2?3l will imply that TV is included in iV^ and that d{Tr,dNr) > ig(-„™Q2,. 
for C (n) large enough. By pulling back the vertices of Nr to T^. (q) under the map exp^ . , we 
can assume that {M'^,g) = (R",eud) (by Lemma [2?3l this operation does not change m and 8 
too much for C(n) large enough). In that case Tr is a real simplex and by convexity argument 
on the distance function, d{Tr,dNr) is bounded from below by the infimum of the distances 
between disjoint faces of a n-simplex of T^-. For a 7i-simplex with mesh m and thinness Q, 
an easy computation, based on multi-linearity of the determinant, gives that this distance is 
bounded from below by (g^„™Q2T. ■ 

We easily infer (3) from what precedes and from the second axiom of the definition of the 
geodesic triangulations. We also have (J U^^Qi?^(A") both closed and open in M, and so 

o-6S„ 

equal to M. This implies that the sets ((1 + il)Ta)^^g cover M. 

The property (4) is obvious when a C r or r C tr and follows from (3) when cr n r = 0. In 
particular, (4) is true when cr or r is a vertex. We now suppose that r and a intersect and no 
one is a subset of the other. As for Point (3), we pull back cr and r in T^M, where z is a vertex 
of (J n r. By Lemma [2.31 it remains to show that (4) is satisfied in the Euclidean case. 

Let a{n) > such that for any fc-face cr and any face t of A" we have 

B((7 \B(crnr,a'=), nB(T,a''-+i) =0 

By dilation based on a vertex of r n cr and rate r < 1 we get 

B{a \ B{a D r, a'^r), a'^+V) n B{t, a'=+V) = 

Since the linear map which maps A" to any for s G S'„ is auto-adjoint with eigenvalues in 
|- C(ra)m ^ C(ti)to] and by Lemma [2.31 we get point (4). □ 

Given a geodesic triangulation T of M, we set, for any x G M, dimT(a;) = inf{|j > 0/3a £ 
Sp,d{x,T^) < This is well defined by Point (2) of Corollary [2H as soon as 

iTT-T < „, NrfL<i A2 • ^i^y simplex cr S S'„, we set 



Sa=B\Ta, — mTj\ IJ B\Tr, ruT 



6m ' ~; ^ 5m 

TCacr 
Sa^\J Sr. 

The following properties follow readily from Lemma 12.31 and Corollary 12.41 It fundamental for 
our application to spectral approximations. It says that, even if a geodesic triangulation is not 
an actual triangulation of the manifold, you can decompose the manifold into some thickening 
of the generalized faces Tg-. 

Proposition 2.5. Let {M'^,g) be a manifold with Sm |cr| < and T be a geodesic triangulation 
of AI with lOmr < inffi ivf, /y , ). We have the followinq properties 

(1) /or any a G S^, we have Sa = {x £ T^j dimT(a:) = "p], 

(2) M is the disjoint union of the {Sa-)a-eK and for any (cr, r) € K^, we have SaC^ Sr = 

SaDT ■ 

(3) for any a e Sp and t e Stn{a) we have that Vo\S„ < C(n) (^)""^ Vol S"^. 
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2.3. Construction of good geodesic triangulations. Given a compact manifold, we can 
use the Riemannian exponential maps to construct some geodesic triangulations at the neigh- 
bourhood of any point with bounded thinness and arbitrarily small mesh (image of some Eu- 
clidean triangulations of the tangent space) and then adapt the Cheeger-Miiller-Schrader ([6]) 
procedure to interpolate these local triangulations in a global, controlled triangulation of the 
manifold. 

To make easier the control of the thinness in our construction, we will work with an alter- 
native (fortunately equivalent in bounded geometry) thinness Qt of triangulations. In that 
purpose we replace the Gramm matrix A"^ by the matrix 



A"" = 



in the definition of the thinness given in the introduction. By Theorem 1 2. 21 if (M", 5) satisfies 
^mW\ ^ and niT < inf (u/, 2a) ^^^^ '^s have 



(i-A2(p:)^)^Ki^,<d(x,.,^(,))^<(i + A^(fi)^)^Ki^^ 

(l-A^(f^)^)^KlL<'^(^'^^^.(0)'<(l + A^(?^)')'KlL' 



(1 - t.\^f)'V, - vl\l^ < d(.,„(,,x,„(,))^ < (1 + A^(5^)^)%r - vl\l^. 



which easily gives A'^i — A^^ ^ 6m'^{Aj^)'^ . This easily implies the existence of some functions 

Ci, C2 such that Qt ^ Ci{<dT,n) and Qt ^ C'2(0T,n-) as soon we have itlt < inf(jM, 
Note that this two thinnesses coincide for Euclidean simplicial complexes and the our thinness 
is essentially the inverse of the fatness used in [5] . 

Let {xi)i(zi be a maximal family of points of M such that the balls Sj;. (lO^/e) are disjoint. 
Let /i,--- , /fe be a partition of / into (non empty) parts such that each (i?^;. (40-ye))ig/^- is 
a maximal family of disjoint balls among the (Bx; (40-y/£))ig/\uj^^^/^ . Since for S Ik and 
for each j < k, Bxi^{4:0y/e) has to intersect at least one ball . . (40-ye) with e Ij, the 
Bishop-Gromov inequality gives us 

^ ^ YolB^. {160 Je) 

fcVolS,,J120V^) < ^VolB,.^(160Vi) ^Vol5x,^.(10Vi)max ^^^^ ' 

and so we have k ^ C(n) for any e ^ 

By iteration, we will construct a family of geodesic triangulations, possibly with boundary, 
Ci, • • • ,Ck in M with mesh less than eC{i,n), thinness ^ C{i,n) and whose vertices of 
the boundary are outside the set ^iGij Bxi{30y/e ~ c{i,n)e) for any e ^ C{iMin, j^). For 
i = k and e ^ C{iM,Ti, j), C'fe will be a geodesic triangulation of M (without boundary since 
Ui?2,. (20v^) = Af) with mesh less than C{k,n)e and thinness less than Q{n) — C{k,n). 

For any e > 0, there exists a constant C(n) > such that K" admits a triangulation with 
mesh less than e and thinness less than C{n). For any i G I, let Ti be such a triangulation of 
Tx^M and T/ the subcomplex whose simplices are those of Ti contained in Bq^, (30-ye) C T^-M. 
For any j £ {1, • • • ,k}, we set Kj the simplicial complex Uig/^.T/. Kj is naturally identified 
with an abstract simplicial complex of M with vertices {exp^.{y),y £ T/,i G Ij}. By Theorem 
[O we have niK^ s$ s/2 and Ok, 2C(n) for any e ^ C{n,iM, 

We set Ci = ^1. Assume that d is constructed. We now construct C^+i by interpolation of 
Ci with i^i+i. For any I G /i+i, we consider in T^^M the complex T/ and the complex Si whose 
vertices are the pull back by exp^^ of the vertices of Ci that are contained in Bq^ (iOy/e), and 
whose simplices have the same combinatorial than in Ci. Using Theorem 12.21 as above we get 
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that Si is an Euclidean complex with mesh less than 2C{i,n)e and thinness less than 2C{i,n) 
for e ^ ^^J-^inf(jMj ^) for r]{i,n) large enough. In what follows, for any Euclidean complex 
T, we denote by T its support, i.e. the union of its simplices. Let Ai (respectively A[) be the 
complex formed by the simplices of T[ (respectively not) contained in Si, and Bi (respectively 
B'l) he the complex formed by the simplices of Si (respectively not) contained in T';. We set 
Di the set of simplices of T^^M that are intersection of a simplex of Ai and of a simplex of 
Bi. Let El be the set of the simplices a of A'^ disjoint from BJ and whose intersection with Di 
is either empty or a union of faces of a. Similarly, we set Fi the set of the simplices a of B'l 
disjoint from Aj and whose intersection with Di is either empty or a union of faces of a. 

Di is a polyhedral complex. We obtain a triangulation D; of Di by first barycentric subdi- 
vision as follows. Let C be a cell of Di. It is a closed, convex polyhedral cell, as an intersection 
cr n cr' with cr in Ai and a' in Bi . For each face a a C 9C, we set pa the isobarycentre <Ja ■ The 
simplices of D'^ are those spanned by all sets Pai, . . . , Pq j , with ctq . C ctq ^ for any 1 ^ z ^ t — 1 . 
By the proof of Lemma 6.3, p. 439-440, and by Lemma 7.1 3) of there exists some con- 
stants /(6, m, n), g{Q, n) and /i(9, n) such that, up to a move of the vertices of T/ by at most 
e2C{i,n)f{2C{i,n),n), the thinness of D'^ is less than g{2C{i,n),n) and its mesh less than 
eh{2C{i, n), n). Actually, we perform this deformation of the complex T/ before the definition 
of the complexes A;, A\, Bi, B[, Di, Ei and Fi. Since 6; allows a control of the thinness of the 
first barycentric subdivision of any intersection cell of a simplex of (the deformation of) T/ and 
of a simplex of Si , there is no circular definition. 

To extend this triangulation to Ei U F; , we keep unchanged the simplices of Ei U Fi that do 
not intercept Dj and we subdivides all the n-simplices of Ei (or Fi) with non-empty intersection 
with Di. For such a simplex a of Ei, we have a r\'D\ — da O dT)\, since any cell of Di is covered 
by some simplex of A[, and so different from a. The triangulation D'^ induces a partition of 
crnDi, which by definition of Ei is a triangulation {da)' of da. We then subdivides a by forming 
all the simplices spanned its barycentre and by a face of {da)'. Thus we get a triangulation Gi 
of El U Fi. Once again, by Lemma 7.1 3) of f5], the thinness of G; is less than g{2C{i,n),n) 
and its mesh less than h{2C{i,n),n). 

We set C^+i ; = GiUD'i and C{i + l,n) = ma.x{l, g{2C{i,n),n),h{2C{i,n),n)). Then C^^^ ^ 
is a simplicial complex of T2;,M with mesh less than eC{i+l, n) and thinness less than C{i+1, n). 
Indeed, since Gi and D[ are simplicial complexes, we just need to remark that a simplex of Gi 
and a simplex of D'^ intersect along a simplex of dD'^ by what precedes. Moreover, any point 
of TJ U Si is either in Di or contained in a rt-simplex of or of i?;'. So (T{ U Si) \ C[^^ j is 
covered by the simplices of AJ (respectively of B'^) that intersect BJ (respectively AJ) or whose 
intersection with Di is not a union of their faces. In the former case, the simplex obviously 
intersects both dSi and dTl and so is at distance from dS\ and from 9Tj less than 2mc'^_^ ^ ■ 
It is the same in the latter case, since if cr is a simplex of (for instance) AJ whose non empty 
intersection with Di is not a union of face of cr, then a face cr' of a must satisfy cr' n Di 7^ and 
cr' \ Di 7^ 0. It first gives that a intersect Si (but is not contained in Si) and so is at distance 
from dSi less than mc'.^^ ^. By definition of Di, a' D Di is covered by some simplices of T/ 
contained in Si, and since T/ is a simplicial complex, we get that a' C Si. So cr' is covered 
by some simplices of Si. Since a' is not in Di, we infer that cr' intersect some simplex of _BJ. 
Hence, a' is at distance from 9TJ less than wc^^ ^ and a is at distance from 9TJ less than 
2mc' ^ J • We infer that (TJ U Si) \ C[_^^ , is covered by the intersection of the 2mc' ^ ^-tubular 
neighbourhoods of 9T/ and dSi . 

We now set Ci+ij the union of the vertices of Ci \ B^i (35v'e) and of the image by exp^^ of the 
vertices of Gi U Dj. We endow it with the abstract structure of complex obtained by gathering 
that of Ci \ Bxi{'iQ\/e) and that of G; U DJ. Since the complex Si is not deformed during the 
previous interpolation and since the only vertices of 5*; that disappears during the interpolation 
are in i?^;, (35-y/e), we really get an abstract structure of simplicial complex on G^+i./ such that 
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any simplex is contained in a n-dimensional simplex. Moreover, by Theorem 12.21 C^+i^; is a 
triangulation with boundary of (M", g) with thinness less than 2C{i + 1, n) and mesh less than 
2C{i + 1, n)e for any e ^ C{n, jm, j)- Finally, by what precedes the vertices of the boundary of 
Ci+i^i are the same as the vertices of the boundary of Ci outside i?^, {30^+C{i, n)e) and are at 
distance less than C(i, n)e from dB^i (30-v/e) inside (30-\/£+C(i, n)e) for any e ^ C{n,iM, j^) 
(once again by Theorem l2.2l) . From this, we get that the vertices of the boundary of Ci+ij are 
outside (u;-^i Upg/^ B^^{30y/e - c{i + l,n)e)) U B^^{30y/e - c{i + l,n)e). 

Ci+ij is just the interpolation between d and T/. But since the family of balls B^^ (40^/e), 
/ G li+i are disjoint, we can interpolate Ci with all the T/ {I S h+i) simultaneously to get 
Ci+i. Note that the constant c{i + 1, n) and C{i + 1, n) will be the same whatever the cardinal 
of /i+i is since the operations done during the interpolation of two different T/ do not interact. 
So we get the geodesic triangulation Ci+i with all the needed estimates. 

Note that the image of any simplex of our geodesic triangulation by the barycentric co- 
ordinates map associated to its vertices gives an embedded simplex of M and thus a true 
triangulation of M whose edges are minimizing geodesic segments. 



3. Estimates on the eigenfunctions 
The following proposition gives bounds on the gradient and Hessian of the eigenfunctions. 

Proposition 3.1. Let {M,g) be a compact Riemannian manifold with delta\j\a\ < A^. For 
any f d Ep, we have that 

(3.1) ll/lloo<C(n)e^^(l + 5l,Ap)t||/||2, 

(3.2) Sm ||d/||oo < C{n)e'^^{l + Sl, Sm \\dfh < C(n)e^A(l + 4, Xp)^\\fh, 

(3.3) SM\\Ddf\\^<Cin)e"'^il + SltXp)'^\\df\\^<Cin)e'""'^{l + Sli\pf^\\df\\,. 

Remark 3.2. The three first inequalities of Provosition \3.1\ are valid under the weaker assump- 
tion (5|jRic > — A^ g. 

Proof. The proof of Proposition [Oi ls based on a Moser iteration schema. For any / G H^{M), 

we have the Sobolev inequality ||/||_22i_ < C{n)e^^ 5m \\df\\2 + II/II2 (if M is a surface, this 

inequality, and what follows, is valid with n — A). 

Let A = D*D be the rough Laplace operator on covariant tensors of (M", g) (where D* is the 
adjoint of D). The rough Laplace operator coincides on the functions with the usual Laplace 

operator and, on the 1-forms it is related to the Hodge Laplace operator by the Bochner formula 

A = A + Ric. For any tensor T on M it satisfies the equality (AT, T) = \DT\^ + ^A{\T\^). 
Let T be any tensor on M and V be a field of symmetric endomorphisms on the covariant 

tensor bu ndle of M . Fix a real V>0 such that <V{T),T> > -V\T\^ for ah tensors T. We 

set u = -\/|Tp + e^, then we have that 

uAu^^A{u^) + \du\^ < iA(|T|2) + \DT\^ = (AT,T) < \(A + V)T\u + Vu^ . 
This inequality and the Green formula gives, for any real fc > 1/2 

,1^. fcM,2 fc' f {du,d{u^''-^)) e f {uAuju^'^-^ 



2k - 1 Vol M 2k -1 Jj^j Vol M 

<^^[m+v)T\u\\u\\f,-'+nu\m 
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We apply the above Sobolev inequality to the function and make then e tends to 0. This 
gives us the following inequality 

(3.4) llTlli^ < (|1T|12\ + ^^'^2\^!!fV ll(^+^)^ll^^ll^ll2r'+a|T||^^) 

We have Ep = Vect{/i, i < p}. For any A: > 1, we set Ak = sup ||'{]]^ (resp. Bk — 

sup H ). Since Ep is stable by A we have 
feEp\{0} Mh 

||A/||2fe < A^kW^fh < A^kXpWfh, W^dfhk < B2k\\Adf\\2 < B2kXp\\df\\2. 

Hence, by applying inequalitv 13.41 to T = f and = (resp. to T = df and V = Ric), for any 
/ G Ep, we get 

— vT=i — ) 

A I. . \ l/fc 



B^<{1+ ^^J^^l ^ ^{n - 1)A^ + 6j, Xp ] B2U 



1-2 



We multiply the inequalities obtained by setting successively k = with v — -^^^ > 1 and 
j G N. Since Am tends to A^o (resp. Em tends to Boo) when m tends to cx), we get 



max(Aoo,-Boo) < n f 1 + - + 6m Xp 

To get a more convenient upper bound, note that 1 + aVb < \/l + b (1 + a) and that the infinite 

1 

product + ^■2yi-i ) converges, hence 

max(Aoo,-Boo) < C"(n)eTA(i + _ i)a2 + Ap)^. 



This gives us the first three inequalities of proposition 13.11 For any / G Ep, we have (see for 
instance [I]) 

(3.5) {KDdJ, Ddf) < {{D*R)df, Ddf) + ^^^^\Ddf\^ + {DAdf, Ddf). 

Om 



Now, if we set u — ^J\DdJY^^ we have uAm < {ADdf, Ddf) . From Lemma 13.51 we infer 

U2 ,r<f„\\2 



(3.6) <^f^^\\u\\ll+ I {DAdf,Ddfy^'-'^+ I {D*Rdf,Ddf)u^^^-'^' 
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The divergence theorem appHed to the field u^^*^ ^^Adf, D,df)'^ , gives (Vfc > 1) 

/ {DAdf,Ddf)u^^''-'^ ^ [ |Ad/|2u2(^-i)-2(fc-l)V(Ad/,M/(eO)dw(e,)«2fc-3 
Jm Jm j 

< / |Ad/|2M2(fe-i) + 2(A: _ 1) / \Kdf\\d% 

JM JM 



.U\U 

M Jm ' 

<^ f |dii|V('=-i) + (2fc-l) / |Ad/|2u2(fc-i) 
2 Jm jm 

<^ f Mu|V('=-i)+2(2fc-l) / (|Ad/|2 + |Ric(d/)nu2(fe-i) 
2 -Im jm 



<^ f |du|V('^-i)+2(2fc-l)(||Ad/||L + ^5^IM/lli 

^ JM om Jm 

The same process applied to vP''^^^'^^ {tvi ^{{R(^, ,^df,D,Ddf)))'^, combined to the equality 



gives: 

{D*Rdf,Ddf)u^^''-^^ 



M 

= f i|M/|V('=-i) + 2(fc-l) V / {R^{e,,e/)df,D,^df)du{e,)u^''-^ 
Jm ^ Jm 



M JM 



< 



2 

fc- 1 



4 

-I 



Jm' ' ^ ' 



Since /^^ |dupM2(fe-i) = ^ J^^ inequality ([311) implies: 



II II r4 

Now, since ||Ad/||oo < Soo||A(i/||2 < -BooAp||d/||oo, we get 



We can now apply the Sobolev inequality to get: 

WDdft^^ < \\Ddm, + kC{n)e^-' \\Ddf\\1-\ (\\Ddf\U + ^2^1(1 + ^l, Ap)i+t) 
< \\Ddf\\lzl2\\Ddf\\^(l + C7(n)e^"'A:(l + (1 + ^v)"^ J^^^^)) 
< \\Ddf\\^^l,\\DdfU(l + C(n)e^"'fc(l + Sl, Xp)^) 
since we can suppose that Sm \\Ddf \\oo > ||<i/||oo- Hence (set k — ai/2 + 1) 

2 

Fi+i< l + kCin)e''"\l + sli\p)^ " Fi, 
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where 



ri-2 



> 1, Fi ^ 



( ) ^'^'^ sequence defined by oq = 2 and 



a;+i = ;^(a; + 2). Since the sequence ai/i/'' tends to n, we get 



WDdfh 

The previous inequality gives 



i=0 



\\Ddf\\oo<C{n)e'^^{l + Sli Xp)^\\Ddf\\2 

But if we integrate the Bochner formula {Adf,df) — iAjd/p + \Ddf\^ + Kic{df,df) we easily 
get 



So we have 



5m WDdfh < ^Jsli Xp + {n-l)A\\dfh 



Sm \\Ddf\\^ < C7(n)e"'^(l + Sti Xp)^\\df\\ 



□ 



Proposition 13.11 implies that at small scale, the eigenfunctions are almost affine. 



Lemma 3.3. Let (AI,g) be a compact n-manifold which satisfies S'^\a\ < A^, and T be a 



geodesic triangulation of M such that lOmx < inf(iM 
function on 10 ■ by 



Sm ' 
2A . 



For any cr G 5„, we define a 



Then we have the following estimates on 10 • Tg- 



11/ - lU < C{n)e^-'Hl + S'm Ap)^ 

OM 



Proof. We set v = X]?=i ^j^J ^'^'^ — ^^Px„{t'^)^ then we have that 



f ilil)) ~ f{X,) ~ g-^D fix), j{0)) \ = |//'^^°^(^)^*^^ 



< 



|i?d/|o7(t)dids |7(0)P < 



This implies that 



i=i j=i 



If we set 9j = Sjk , then this inequality gives that 



<Cin)e-^'\l + 6lX,)'^{^r\\f\U. 



which combined with the previous inequality gives the first result. 

Set L the linear form on Tx^M such that LI = L o exp^^^ . For any w € Tx^M we set w{t) 
the parallel transport of w along 7. Then by Theorem 12. 11 we have that 



(3.8) 



dLt{w{l)) ~ L{w)\ < \\L\\\dexp-^l{w{l))-w\<i\\L\\\w\{^r. 



14 



ERWANN AUBRY 



On the other hand, we have that 



\df{wil))-dfiw)\< I -df{w{t))dt 



< 



Ddf{j,w{t))dt < \\Ddf\Uv\\w\. 



Hence we get 



(3.9) \{dLl ^ df)w{l)\ < (2\df-L\ + C{n)e-'''{l + S^\)'^\\df\\^j)\w\. 



Now, by InequaUty (|3.7I) we have that 



and so 

\df - = - L){vJ)idf - < \\A-^ ^[(d/ - L)(v^)r 

jk j 

<enn)e'-'^ii+si,\r\'^ndf\\i. 



If we combine this inequaUty with Inequahty (I3.9p , we get 

\dLi-df\<Cin)e'e'^'^{l + S'Xp)'^j\\df\U 
which gives the result by Proposition 13. II □ 

4. Proof of Theorem 11.11 
4.1. A discrete norm. We prove that gives an approximation of II/II2 on Ep. 

Proposition 4.1. Let (M'^,g) be a compact Riemannian manifold with S%j \cr\ < and T a 



geodesic triangulation with lOmr < mi(iM, 



5 m 



C(n)Ae? 



Then, for all f £ Ep, we have 



M 



5m ' 



A/ 



Proof. Corollary HjH (2) implies that 

/ [Lif + (2II/II00 + \\Li- f\\oo)\\Li - f\U By LemmaESland 

J(l+v)T„ 



We have 



Proposition 13. 1[ we have 



(l+'7)T„ 



{Lif + C{n)e'-'\l + 5l, Xpf^' (?)^||./||^ Vol((l + r?) • T.) 



dv, 



By Theorem 12.21 we have that 

[f[X,) + - fix,)]) dv, 



< (1 + 4A^(-)^)' 





(i+';)A„ 



< 



< 



2VdetA 



f{X„) + Y,e,[f{x,^^,))-f{X^)\) dz;,,^ +2"4||/||L(A'(y)' + r;)VolA, 

Am, 



^— ^ E + C(n)e2"( — Vol(l 



0<j<j'<n 
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So, we have 



+C(n)e^"^'^(l + 5l, Ap)'"'e2"(^)2||/||2 Vol((l - ,y) • T.) 

By summing on a e 5„, we get - |i?(/)|2 < C{n)e'^'^{l + Sl, Xpf"'( f)^Q^^ J^^ p. By 

the same way, we get the reverse inequahty. □ 

4.2. A discrete Dirichlet energy. We prove that <7t(-R(/)) approximates \\df\\2 on Ep. 

Proposition 4.2. Let (M",g) be a compact Riemannian manifold with Sm |cr| < and T a 
geodesic triangulation of M with IQmT < inf (z^f , c(^n)e"A ) ■ '^^^''^ /'^'^ '^'^2/ / ^ ^p; '"'^ have 
that 

Jm om Jm 

Proof. As in the proof of Proposition 231 we have 

M/P< / \dLl\' + C{n)e^e'-'^{l + sliXX^^\\df\\lyo\{{l + 7^)-T,). 



Using Lemma ISJ and the fact that J^^ \L\^ = ^^i^ Efc.i=i where L is the 

linear map defined in the proof of Proposition 14. II (i.e. L{v^) = f{xi^(k)) — fi^a)) we get the 
result. □ 

4.3. A Withney map. Wc construct an extending (i.e. Withney type) map which to each 
{Ui) £ associates a function / : A/ — > R such that f{xi) = yi for all 1 < z < iV. This 
function / has to be such that J f^ and / be close to \yi\'^ and qxiyi). In that purpose 
we take / almost linear by part. 

We need first some controlled partitions of the unity on AI associated to the geodesic trian- 
gulations. 

Lemma 4.3. Let {M^^,g) he a compact, Riemannian n-manifold with 5{MY\a\ < A^ and T 
be a geodesic triangulation of M with IOtut < {iM, )A^o<t"^ ) ' There exists some smooth 
functions {(Pa-)aes„ such that 

(1) (fia- ■ M — > [0, 1], (fia- ~ I on So- and Supp ipa C 5*0- for any a G Sn, 

c 2 

(2) E<Tes„ ¥'(7 = 1 and \dip^\ < '^^"•'^| for any a e Sn- 

Proof We set C = ^-TJ^- Since the distance in {Tx,M,gx,) between (1 + C)A^ and Tx,M\ 
(1 + 2^) Ao- is larger than ^^^^ (see the proof of Corollarv l2.4l) and since expj^^ is a 1 + K'^{^Y- 
Lipschitzian map, there exists a function ipa ■ M ~> [0, 1] such that tAo- = 1 on {l + ()Tcr, tpa — 
outside (1 + 2C)To- and \dipa-\ < . We set ipa- — ^ ''^''^ ^ , which is well defined since by 

CoroUarv 12.41 we have Ures„(l + O^t ~ M . By the same kind of arguments as in the proof of 
CoroUaryO we have (1 - lc,)T„ 0(1 + 2C,)Tr = for t^cf, and so (/j^ = = 1 on (1 - 4C)r^. 
We have obviously ipa — ^ outside (1 + 2C,)Trj- By a volume argument, the number of non zero 
term in the sum dxtl'T is bounded from above by 2"92", and so we have that 



< 



5ri2"e 



□ 
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Definition 4.4. Let fx = (vi) G M^. For any a e K , we set L'^ : {9k) ^ F ^J2k ^kUi^ik) ■ For 
rriT small enough, the function L^J' = L'^ o [Bq)~^ is defined on So- and VF(/t) = faLl'^ 

creS„ 

is well defined on M . 



This extending map W satisfies the foUowing properties. 

reposition 4.5. If lOmr < M{im, 3^;;^X(I?A)ep") ^^^'^ 

(1) RoW = IdRN, 

(2) for any fx G K^, we have that 

\f m/rf -|(/t)||| <C(n)(l + A)e2"(^)2(|/^|| + ^2^qj,(/^)), 
Jm 

(3) for any (fx) G M.^ , we have that 
\dWifT)\^ - gT(/T)| < C(n)e20"^^<ZT(/T) 



\dW[fT)r -qTUT)\ < C'(njH-""'"' 

M 

2 



Proof. Point (1) is obvious. We now prove point (3). By Inequality p.8p . we have ||<ii^^ 
\L.\^\ < 16|L,|2(^)'. Since 

n 

\FA^ = XI i^<yf^iyiAk) - yv(o))(y»„(/) - yv(o)), 
fc,;=i 

we have 

> X VolT.(l - 16(^)^)"+^(1 - ^r\L.\' > (1 - C(n)^),H/^). 
On the other hand, we have that 

n — l 

n—l p 

<{l + C{n)j)qT{fT) + Y.T. \dW{fT)f, 

and for any t £ K 

[ \dw{fT)\'^ [ I E ^(i^"'/'-)!' 

<2[ { J2 V.\dLt-\T + '2 [ I 5] ^^-d'/'.r. 
Since < (1 + 8(^)2)|L^|, and by Proposition O we have 

n—l „ n—l . 

P=OtGSp •'^^ creSt„{T) p=QTeSp <yeSt„{T) 

< C{n)^ V |L.pVol5. < C{n)'^qT{fT). 
OM ~rz OM 
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To bound the remaining term IX^o-est (t) ^cr^ '^'fi'^l^ ' "^o G Stn{T) and tq the reah- 

sation of t associated to expj^^ . Then there for any x Sr, there exists x' G tq such that 

d{x,x') < ^aP+\ Hence we have |L;^^ (x) - (x')| < 2|L^J2™!q,p+i. Note that the barycen- 
tric coordinates {9k) of x' in the simplex a satisfy 0^ = if Xi^ (^j ^ r. For any a £ S't„(T), 
we set Xa the point whose barycentric coordinates in a are 9' , , = 9 -i, > if x is a vertex of 

T and = otherwise. By applying Lemma [2.31 at most twice, we get d{x',Xa) < 200A^^. 
Since Lfj'{xa) — i^^(x'), we have that 

n-ix) - Lil{x)\ < \Lf/ix) - Ll-ix,)] + \Lllix') - Lllix)\ < + 
Since J^a '^Va — 0; Lemma 331 gives us 

n— 1 « n— 1 « 

EE / I E i^^^^^r = EE / I E iLff - L^JM^A" 

n-1 

^ E E CWei""'#5t„(r)( 5] |L.p + \L,,\')Yo\Sr 

p=lT6Sp o-GSt„(T) 

We now prove point (2). As in the proof of point (3), we have that 
and 

n—l „ n—1 . 

EE / m/T)p = EE / I E ^^^^1^ 

^EE E C(n)e^"/ |L^-p = C(n)e^"^ /_ |L/-r 



Hence we have that 



and since by Theorem 12. 2[ we have 

(i-(^)^)"i/Tii = (i-(^n"E/^ i^n' 



2 



it only remains to bound from above I^^^P- If we set Q = n9^"a""''^^, then Sa\T„ C 

(1 + C)A. \(l-C)Aa andso 

JS,\T^ J(1+C)A„\(1-C)A„ 
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Let H\ be the dilation in F of factor A and centred at the iso-barycentre of A". By the 
fundamental theorem of the calculus we have 



1 



(!-(" + 1)0" Ai-OA^ J A. 

We set h^{n + l)^C,. Since 

|(L-)2(iJi_(„+i)^(y)) - (i-)2(y)| < 2|i-(y)|||dL^|/^&^+||dL^P4//i' 

<C{n){^f{\L%y)\^ + 5lA\dL^') 

OM 

then we have 

1 



(i^- {Ln' <C{n){ — Y (i-) 



By the same way we get 



jm 



1+C)A„ JA„ OM Ja^ 



Hence E.s„ /(i+c)AA(i-C)a„ 1^1' ^ C^(")(i^)'(l/T|| + <Sif 9t(/t)), which combined with the 
previous inequalities gives the result. □ 

4.4. Conclusion. By Propositions 14 . 1 [ 14 . 2 1 and 1 1 . 71 we can almost bound from below the first 
eigenvalues of (M",g) by the eigenvalues of qr with respect to (•, ■)t- To have an error bound 
that depends on p and not on Ap, we use the following rough version of a well known result due 
to S. Cheng. 

Lemma 4.6. Let i^Al^^g) he a couipact RieTnauuidTi Tnanifold such that I < A^. For any 

k e N, we have Sm ^fe < C{n){^)'^e'^ k'^ . 

We infer the following theorem. 
Theorem 4.7. Let e g]0, 1[ be a real number, {M",g) be a compact Riemannian manifold such 
that Sm |cr| < A^ and T be a geodesic triangulation of (M", g) such that — ^( 



C{n) iM 

then we have Xp{qT) < Xp{M" , g)(l 



3n 



Once again, by Propositions 14. 51 and II. 71 we can bound from above the eigenvalues of {AF\g) 
by the eigenvalues of qx with respect to (•, ■)t- Note that to bound Xp{qT), we just have to use 
Theorem H?71 

Theorem 4.8. Let e g]0, 1[ be a real number, {M^,g) be a compact Riemannian manifold such 

I C(n)i \ 3"^ 

that Sm I'^I < A^ and T be a geodesic triangulation of (M", g) such that < e( j 
Then, we have that 

Ap(M",5) < Ap(<7T)(l + e). 



Proof. When applying Proposition ! 1 . 71 with Ei the space spanned by the p+1 first eigenfunction 
ofgT, Propositionl431 gives a > 1-C{n)e^e^'^if^y {l+Sli Xpiqr)) and /3 < l+CHO^""' j^.. 

Now, by Theorem O and Lemma Ol we have Sli ^p(9t) < '^^"^^f □ 
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5. Proof of Theorem 11.41 



5.1. Approximation of the eigenfunctions. To get the relations between the eigenfunc- 
tions of {M",g) and the discrete eigenfunctions, we first prove the fohowing resuh, where the 
notations are the same as in the introduction. 



C(v n A ^M-) 1 

For anyfe Ep, we have ||i?(/) - Pp o i?(/)||2, ^ y^' '° ij^)^ \\Rif)\\T where Pp is 



Lemma 5.1. Let 5 > Q and assume that Xp{M) + rj ^ Ap+i(M). 

(/: 

the orthogonal projection from M.^ to the space Fp spanned by the first p eigenfunctions of qt ■ 

For any (y,) € Fp, we have \\Wiy,) - Qp o Wiy^^^ < C(p,n, A, ^,r,)(2^)^^|lW^(y,)ll?^, 
where Qp is the orthogonal projection from L'^{M) to Ep. 

Proof. We use the same idea as in 8 . We consider in A^+^M^ the operator A{vo A ■ ■ ■ Avp) — 
SiLo "0 A • • • A Axivi) A ■ ■ ■ A Vp, where is the symmetric operator such that qrix) = 
{At{x),x)t. The eigenvalues of A are exactly the sum Xi-^{T) + • • • + Xi^_^-^{T) with ^ ii < 

■ ■ ■ < ip+i N, and so its first eigenvalue is Ao(r)H hAp(r). We set R : AP+^Ep Ap+V 

defined by R{fo A ■ ■ ■ A fp) = R{fo) A • • • A R{fp). Since |i?(/o A • • • A = (det G)^, where 

G is the Gramm matrix of the family i?(/o), • • • ,R{fp), and since there exists Cp such that 
l(detG)^ ^ 1| ^ Gp|G — /p+i| for G near /p+i, Proposition 14. 1 1 and Lemma l46l give us 

|l^(/o A • • • A fp)\l - 1| ^ G(p, n, e, A, ^-^)C^f 

iM OM 

By Propositions 14.11 and 14.21 and by Theorem 11.11 we get 

p 

\{AR{fo A • • • A /p), i?(/o A • • • A fp)) - ^ A,(r)|i?(/o A • • • A /p)|'| 

P 

= lE(^(/o) ^ • • • ^ - MT)R{f^)) A • • • A Rifp),Rifo) A • • • A i?(/p))| 

p 

= I E (^(/o) ^ • • • ^ (Mm)) - \{T)R{f,)),R{f,))G''R{f,) A • • • A Rifp),Rifo) A • • • A i?(/p))| 

hj=0 

= |f]G"(detG)2(g^(i?(/,)) - A,;(T)|i?(/,)||)| 

4=0 

f]G"(detG)2(gj,(i?(/,)) - \\df,f )\ + |f]G"(detG)2(A,(M) - A,(r))| 



i=0 i=0 

if] G"(det G)2A.(r)(l - < Cip, n, 8, A, 



Let (j/i) an orthonormal family of eigenfunctions of qx associated to the eigenvalue (Ai(T)). We 
set R{fo A ■ ■ ■ A fp) ~ ago A ■ ■ ■ A yp + n with n orthogonal to yo A ■ ■ ■ Ayp. The above estimates 
give us + jnp - 1| <^ C{p,n,e, A, f^)C-f^)^ and S\n\'^ C{p,n,Q, A, l^){f^)^ , from 
which we easily get 

P C(p n <d A ^" ) ( ™^ ) 

mh) - PpiRimir < l^(/o A • • • A /p) - yo A • • • A 2/p|^ — — ' ' ' 



The proof of the other estimate is exactly the same, but we first have to bound from below 
the gap Ap+i(r) — Ap(r) using the bound on the gap Ap+i(Af) — Ap(A/) and Theorem ll.il □ 
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We easily infer Theorem 11.41 from the previous Lemma. Indeed, by Proposition 14.11 the 
quadratic form \Pp o i?p on Eg takes only values less than 1 + C(Jj^)'^ on the unit sphere of 

Eq and so its trace with respect to (•, ■)t is less than p + C{y^)'^ (complete an orthonormal 
basis of the kernel of PpoR). But the previous lemma, applied for the spectral gap at p, implies 
that the trace of \Pp o i?p on Ep is close to p and so C{j^)e^ + J2p+i^i^q \Pp ° ^(/OP — 

tr\PpoR\^ ^p + C{f^)^, and so Ep+i^.^, \PpoR{m^ i'c(f^)i^. This gives the resuh 
when combined with the previous lemma applied to the spectral gap at q. 
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